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In relativistic heavy-ion collisions, anisotropic collective flow is driven, event by event, by the initial 
eccentricity of the matter created in the nuclear overlap zone. Interpretation of the anisotropic flow 
data thus requires a detailed understanding of the effective initial source eccentricity of the event 
sample. In this paper, we investigate various ways of defining this effective eccentricity using the 
Monte Carlo Glauber (MCG) approach. In particular, we examine the participant eccentricity, 
which quantifies the eccentricity of the initial source shape by the major axes of the ellipse formed 
by the interaction points of the participating nucleons. We show that reasonable variation of the 
density parameters in the Glauber calculation, as well as variations in how matter production is 
modeled, do not significantly modify the already established behavior of the participant eccentricity 
as a function of collision centrality. Focusing on event-by-event fluctuations and correlations of the 
distributions of participating nucleons we demonstrate that, depending on the achieved event-plane 
resolution, fluctuations in the elliptic flow magnitude V2 lead to most measurements being sensitive 
to the root-mean-square, rather than the mean of the V2 distribution. Neglecting correlations 
among participants, we derive analytical expressions for the participant eccentricity cumulants as 
a function of the number of participating nucleons, A^part, keeping non-negligible contributions up 
to O (l/A^part)- We find that the derived expressions yield the same results as obtained from mixed- 
event MCG calculations which remove the correlations stemming from the nuclear collision process. 
Most importantly, we conclude from the comparison with MCG calculations that the fourth order 
participant eccentricity cumulant does not approach the spatial anisotropy obtained assuming a 
smooth nuclear matter distribution. In particular, for the Cu + Cu system, these quantities deviate 
from each other by almost a factor of two over a wide range in centrality. This deviation reflects 
the essential role of participant spatial correlations in the interaction of two nuclei. 



I. INTRODUCTION 

One of the strongest pieces of evidence for the for- 
mation of a thermalized dense state of unconventional 
strongly interacting matter in ultra-relativistic nucleus- 
nucleus collisions at the Relativistic Heavy Ion Col- 
lider (RHIC) P-Q stems from the strong anisotropic col- 
lective flow measured in non-central collision events 0- 
[T2j . Studies of the final charged particle momentum dis- 
tributions have revealed strong collective effects in the 
form of anisotropics in the azimuthal distribution trans- 
verse to the direction of the colliding nuclei, and theory 
holds that their anisotropy around the beam axis in non- 
central collisions is established during the earliest stages 
of the evolution of the collision fireball [l3l - fl6| . The main 
component of this anisotropy is called "elliptic flow" and 
its strength is commonly quantified by the second coeffi- 



cient, V2-, in the Fourier decomposition of the azimuthal 
momentum distribution of observed particles relative to 
the reaction plane [l3| . 

By now there exists an extensive data set of elliptic 
fiow measurements in Au + Au collisions at RHIC as 
a function of center-of-mass energy, centrality, pseudo- 
rapidity and transverse momentum [sl-fT^. The mag- 
nitude of the observed flow anisotropy is found to be 
strongly correlated with the anisotropic shape of the ini- 
tial nuclear overlap region. This is expected if inter- 
actions among the initially produced particles are very 
strong, leading to anisotropic pressure gradients, which 
transform the initial spatial eccentricity into a final mo- 
mentum anisotropy [l8| . 

Quantitatively, the connection between initial spatial 
and final momentum anisotropy is explored by hydrody- 
namical calculations that, for a given equation of state, 
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relate a given initial source distribution to the final mo- 
mentum distribution of the produced particles. For Au-I- 
Au collisions at the top RHIC energy, yj — 200 GeV, 
such calculations are in good agreement with the elliptic 
flow data at mid-rapidity 0, [l^. From similar stud- 
ies, it has been numerically established that, for not too 
large impact parameters, the final magnitude of the el- 
liptic fiow is proportional to the initial eccentricity, e, 
used to characterize the spatial anisotropy in the trans- 
verse plane of the matter created in the overlap region 
of the colliding nuclei [H, [lOl . More generally, one ex- 
pects the ratio of elliptic flow and eccentricity, v^jt^ at 
mid-rapidity to be a universal function of density and 
size of the system at the time when the elliptic flow de- 
velops ("f2/e scaling") [2Tl - [23j . In hydrodynamics, this 
function depends parametrically on the speed of sound 
in the fireball medium [20j . 

The elliptic flow in Cu -I- Cu collisions at RHIC was 
found to be comparatively large, especially for near- 
central collisions, reaching almost the same magnitude 
as in Au -I- Au collisions for the same fractional cross 
section [l^, HH > and much larger than expected from hy- 
drodynamical models [2^. A quantitatively meaningful 
comparison of the elliptic flow values measured in Cu-|-Cu 
and Au + Au collisions requires dividing out the differ- 
ence in the eccentricity of the nuclear overlap zone since, 
for a given centrality, the average eccentricity depends on 
the size of the colliding nuclei. For the same size of the 
overlap zone similar densities are achieved in the two col- 
lision systems [13, H^, but the Cu-|-Cu system exhibits a 
significantly smaller spatial eccentricity. If one scales the 
measured by this eccentricity, using its conventional 
definition in terms of the spatial deformation of the av- 
erage transverse distribution of participating nucleons at 
a given impact parameter, one is led to the paradoxi- 
cal finding that the smaller Cu -I- Cu system translates 
the initial spatial deformation more efficiently into a fi- 
nal momentum anisotropy than the larger Au -f Au sys- 
tem mil. 

However, this conclusion depends on the definition of 
the eccentricity used in the scaling procedure. In Ref. [2^ 
it was pointed out that the shape of the nuclear overlap 
fluctuates from event to event, and in Ref. [soj it is em- 
phasized that the orientation of its major axes relative 
to the reaction plane (defined by the beam axis and im- 
pact parameter vector) fluctuates as well. This is illus- 
trated in Fig. [T] For a given event, the actual distribu- 
tion of the participant interaction points in the overlap 
zone can cause the overlap zone to be tilted with respect 
to the reaction plane. The participant eccentricity (epart, 
Eq. (jA6|) [lOl) takes this into account by using the princi- 
pal axes of the overlap zone. Because of the fluctuations, 
the ensemble average of the event-wise participant eccen- 
tricity is not identical with the standard eccentricity (eg, 
Eq. (|A2[) ) of the smooth overlap distribution which is ob- 
tained by averaging the participant density in the over- 
lap region with respect to the reaction plane over many 
events. 




FIG. 1: Schematic illustration of a nucleus-nucleus collision 
depicted in the transverse plane. The principal axes ix' and 
2/') of the area formed by the participants are tilted with re- 
spect to the reaction plane given by the axes (x and y) of the 
transverse plane. 

Since hydrodynamic collective flow is not a property of 
the event ensemble, but rather develops independently in 
each collision event, its driving force is the shape and de- 
formation of the initial distribution of produced matter 
in each event. To investigate the validity of the hydro- 
dynamically predicted V2/e scaling one should therefore, 
in principle, compute the ratio W2/epart event by event, 
before taking its ensemble average (w2/epart)- This is, 
unfortunately, not possible in practise since the initial 
spatial eccentricity tpart of a given collision event cannot 
be measured, and a statistically accurate determination 
of the elliptic flow V2 also requires summing the hadron 
momentum spectra over many events, so only its ensem- 
ble average (^2) is known. In practise, the best way to 
approximate (w2/epart) is to scale the measured elliptic 
flow {V2) by a calculated average participant eccentricity 
(fpart) or by a higher moment of it (see below). 

In Ref. [2^ . [3OI the participant eccentricity scaling was 
studied using Monte Carlo Glauber (MCG) calculations, 
where Cpart is computed for each event from the trans- 
verse position distribution of nucleons participating in 
the collision, taken in its individual major axis frame. 
For large nuclei, event- wise fluctuations in the transverse 
density distributions are small, except for the most pe- 
ripheral collisions. Nonetheless, as one approaches zero 
impact parameter (i.e. in almost central collisions where 
both eg and epart are tiny), these small density fluctua- 
tions still cause significant relative fluctuations of epart, 
resulting in a non-negligible difference between the par- 
ticipant and the standard eccentricity, even in Au + Au 
collisions. For the smaller Cu + Cu system the event- 
wise fluctuation effects are much stronger and seriously 
affect the eccentricity over the entire range of impact 
parameters [13]. The participant-eccentricity-scaled el- 
liptic flow {V2) I (cpart) thus differs appreciably from the 
standard eccentricity-scaled elliptic ffow. It appears that 
scaling V2 with the participant eccentricity unifies the 
eccentricity-scaled elliptic flow {v2j / (cpart) across the 



3 



Cu + Cu and Au + Au collision systems [23, HOl, even 
differentially as a function of transverse momentum and 
pseudo-rapidity [3l[. Furthermore, first measurements of 
elliptic flow fluctuations have recently been reported in 
Au + Au collisions at ^^s^^ = 200 GeV [H, ^ . The 
relative fluctuation magnitude o'ep,i,t/{^pa-rt) from MCG 
is in striking agreement with crv2/{'^2) from data [3^ . as 
expected if initial state fluctuations are combined with 
hydrodynamic evolution. 

The initial success of the participant eccentricity calcu- 
lated in the MCG approach immediately suggests a new 
set of questions: 

• How robust are the participant eccentricity results 
to the parameters characterizing the nuclear den- 
sity distribution (radius, skin depth and nucleon- 
nucleon potential)? 

• What is the effect of varying the assumptions about 
matter production (locality, participant and binary 
collision weighting)? 

• What features of the MCG initial state distinguish 
it from the usual optical Glauber model picture? 

• More specifically, what is the impact of the fluc- 
tuating initial conditions on the suggested [H, [2§| 
use of cumulant approaches? In particular, which 
moment of an underlying fluctuating flow distribu- 
tion is measured by the (standard) event-plane flow 
method [l3? 

These questions will be addressed in the present pa- 
per. In addition, following and improving on Ref. [23| . 
we derive analytical expressions for the eccentricity cu- 
mulants in terms of moments of the initial spatial matter 
distribution, including all leading terms. Furthermore, 
by comparing with the numerical MCG model, we show 
that the analytical expressions are misleading as they ne- 
glect important effects arising from spatial correlations 
between the participating nuclcons. 

II. MONTE CARLO GLAUBER MODEL 

To estimate the geometrical configurations of colliding 
nuclei, one typically constructs models based on rather 
generic assumptions about the constituent makeup of a 
typical nucleus. In this context, it is fairly standard to 
assume that nuclear matter in a nucleus is distributed ac- 
cording to the charge distributions seen in electron scat- 
tering experiments. There are two ways of expressing 
these densities in actual calculations (for an overview, see 
Ref. [s^ l and references therein). One way is to assume 
a smooth matter density, typically described by a Fermi 
distribution in the radial direction and uniform over solid 
angle (in the case of spherical or near-spherical nuclei), as 
done in "optical" Glauber calculations [s^, [s^ . It should 
be noted that this method neglects some potentially im- 
portant correlations between participating nucleon posi- 
tions as will be discussed further in Section [IV CI 



A related, but fundamentally different approach is to 
distribute, event-by-event in a stochastic manner, indi- 
vidual nucleons according to the smooth matter distribu- 
tion and to evaluate the collision properties of the collid- 
ing nuclei by averaging over multiple events using Monte 
Carlo methods [33, [3g|. The key ingredients in MCG 
calculations arc the following: 

1. The nucleon position centers in each nucleus are 
distributed according to a probability distribution 
given by the smooth nuclear density function, p. 
One can think of the smooth nuclear density as a 
quantum mechanical single-particle probability dis- 
tribution for the nucleon positions and their actual 
values in an individual collision event as a "mea- 
surement" of their positions in a given collision ex- 
periment. 

2. The nucleons are assumed to travel in straight tra- 
jectories along the beam direction throughout the 
reaction, i.e. their transverse positions are "frozen" 
during the short time when the two high-energy 
nuclei pass through each other. 

3. The nucleons interact with nucleons in the oncom- 
ing nucleus by means of the nucleon-nuclcon in- 
elastic cross section (ctnn) appropriate for the beam 
energies under consideration (measured in proton- 
proton collisions). The nucleon-nucleon collisions 
occur and produce particles independently, i.e. dy- 
namical correlations among the nucleon positions 
in the multi-particle nuclear wave function are as- 
sumed to be negligible. The only correlations in the 
model are of geometrical nature and due to clus- 
tering effects from the interaction process itself, as 
explained in Section HV CI 

Commonly, a nucleon-nucleon collision in the reaction 
is defined to occur if the Euclidean transverse distance 
between the centers of any two nucleons is less than the 
"ball diameter" , 

D = ^o-nn/tt. (1) 

More specifically, the steps of the PHOBOS Monte 
Carlo 3 calculation for a single event are the following: 

• Impact parameter selection: The two nuclei are 
separated in the x-direction by an impact param- 
eter, b, chosen randomly according to dN/db cx b 
up to some large maximum (femax — 20fm> 2Ra)- 
Thus, nucleus A is defined to be centered at 
{x,y} ~ {—6/2,0} in the transverse plane, while 
nucleus B is centered at {x,y} = {+6/2,0}. In 
addition, both nuclei are centered at z = 0, since 
the longitudinal coordinate of each nucleus does not 
matter for the subsequent steps ^. 



^ Throughout the paper, we will keep the common choice that the 
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• Makeup of nuclei: For each nucleus, wc loop over 
the number of nucleons, Nj^ and Ng, and for each 
nucleon center point choose random, uniformly dis- 
tributed azimuthal and polar angles, as well as a 
radius sampled randomly according to the radial 
density distribution p{r). Additionally, to mimic 
excluded volume effects, one may require a mini- 
mum inter- nucleon separation distance (dmin) be- 
tween the nucleon centers of all nucleons in the 
nucleus. This introduces a geometrical correlation 
among the nucleon positions. In the construction 
of the nuclei, we make sure that the center-of-mass 
of the nuclei is correctly positioned, i.e. we achieve 
Y^Xi = ±6/2, ^j/j = and X) = by shifting 
all nucleon centers by the average offset determined 
after the positions of all nucleons in the nucleus 
have been generated. Thereby we ensure that the 
nuclear reaction and the otherwise arbitrary MC 
frames coincide. 

• Collision determination: For each nucleon in 
nucleus A, we loop over the nucleons in nu- 
cl eus B. If t he 2-dimensional Euclidean distance 
\/ Ax^ + between the nucleon from A and the 
nucleon from B is less than D as defined in Eq. ((T|) , 
the number of collisions suffered by both nucleons 
is incremented by one. If no such nucleon-nucleon 
collision is registered for any pair of nucleons, then 
no nucleus-nucleus collision occurred. Counters for 
determination of the total (geometric) cross section 
arc updated accordingly. 

Having access to the number of collisions suffered by 
each nucleon according to this purely geometrical (clas- 
sical) prescription allows straightforward calculation of 
-/Vpart, the number of nucleons which are struck at least 
once, and A^coii, the total number of nucleon-nucleon col- 
lisions. The latter is defined as the sum of collisions suf- 
fered by nucleons in one nucleus with nucleons from the 
other one (to avoid double counting). For every collision, 
the calculation keeps track of the position and status of 
each nucleon in the event, for later usage in the calcu- 
lation of the spatial eccentricity or any other interesting 
quantity. 

The default parameters of the PHOBOS Glauber cal- 
culations for Au -|- Au and Cu 4- Cu collisions at ^Snn = 
200 GeV are listed in Table [D The nucleon-nucleon in- 
elastic cross section of a^N = 42 mb is from Ref . [1^ , 
while the parameters for the Fermi distribution, 

p{r) cx ( 1 + exp (- — -]] 



reaction plane, defined by the impact parameter and the beam 
direction, is given by the x- and z-axes, while the transverse 
plane is given by the x- and j/-axes. 



TABLE I: Default parameters used in the PHOBOS MCG 
calculations for Au + Au and Cu + Cu at v/Snn ~ 200 GcV. 



System 


(Tnn [mb] 


Na/b 


i?[fm] 


a [fm] 


dmin [fm] 


Au + Au 


42 


197 


6.38 


0.535 





Cu-I-Cu 


42 


63 


4.20 


0.596 






TABLE II: Baseline, minimum and maximum values of pa- 
rameters in the PHOBOS MCG calculations for Au-|- Au and 
Cu-|-Cu collisions used in sensitivity studies of the eccentricity 
definitions. 







Au + Au 






Cu-f Cu 




Parameter 


Base 


Min 


Max 


Base 


Min 


Max 


(Jnn [mb] 


42 


30 


45 


42 


30 


45 


R [fm] 


6.38 


6.25 


6.51 


4.22 


4.14 


4.30 


a [fm] 


0.535 


0.482 


0.586 


0.596 


0.536 


0.656 


dmin [fm] 


0.4 





0.8 


0.4 





0.8 



i.e. the nuclear radius R and the skin depth a, are from 
Ref. [40j . The minimum inter-nucleon separation dis- 
tance is set to dmin = fm, i.e. we generally ignore geo- 
metrical correlations in the multi-nucleon wave function 
due to a hard core, since their effect, especially on the 
participant eccentricity, is found to be small, as we will 
report in Section fill Al 

A few comments are in order explaining what this 
model delivers and how we subsequently interpret its 
output. As described so far, the MCG model records 
only the (transverse) position and collision status of each 
nucleon. No particles are produced in the calculation 
and dynamical correlations among the nucleons in the 
nuclear wave function are neglected. By specifying the 
nuclear positions exactly, i.e. as long as we do not allow 
for a smearing around the points given by the model, 
we are prohibited (by quantum mechanical uncertainty) 
from imposing any constraints on the momenta of the 
scattered nucleons and, in consequence, of the particles 
produced by the collision. Source eccentricities calcu- 
lated directly from the distribution of (exact) nucleon 
positions obtained from the MCG model can therefore, at 
least in principle, not immediately be assumed to repre- 
sent the eccentricity of the produced matter distribution 
which drives the anisotropy of the subsequent collective 
expansion. For this reason, we extend in Section IIII Bl 
the model by smearing the resulting nucleon-nucleon col- 
lision points with a profile function in order to model 
the production of (approximately thermalized) matter 
with finite temperature and restricted particle momenta 
in the neighborhood of the collision points delivered by 
the MCG model. The (in-)sensitivity of the source eccen- 
tricity of these matter distributions to the parameters of 
the smearing profile is studied in detail. 
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Number of participants Number of participants 




Number of participants Number of participants 

FIG. 2: Participant (top) and reaction plane (bottom) eccentricity in Au + Au (left) and Cu + Cu (right) collisions as a function 
of ATpart for variations of the Glauber parameters as given in Table |IT1 The inset shows the ratio with respect to the baseline 
calculation in each case. 



III. ROBUSTNESS OF THE ECCENTRICITY 

In this section, we evaluate the effects of variations 
in the nuclear density distributions and of various as- 
sumptions about the sources and spatial localization of 
the initial matter distributions on the eccentricity and 
its centrality dependence. The two definitions of eccen- 
tricity considered in this section are the reaction plane 
eccentricity (see Eq. (jA3[) ) 



and the participant eccentricity (see Eq. (|A6[) ) 



y2 _ ,2^2 



Cpart 



where ct^, ct^ and Gxy are the (co-) variances of the 
participant-weighted nucleon distribution in a given 
MCG event ^ . Their definitions and relation to the stan- 



Both definitions have already been used in Refs. I24l 13' 



dard eccentricity (eg, Eq. (|A2[) ) of the event-averaged dis- 
tribution, are discussed in Appendix [X] 



A. Variation of Density Parameters 

Before Refs. [H, [13 the purpose of MCG calcula- 
tions was to estimate global properties of nucleus-nucleus 
collisions, i.e. to calculate centrality- and eccentricity- 
related quantities on average^ based on (large) samples 
of Glauber events. Since both the participant eccentricity 
and the reaction plane eccentricity, explicitly involve an 
interpretation of each MCG event individually, it is im- 
portant to understand their dependence on the choice of 
the MCG calculation parameters. A number of sources of 
systematic error are studied by varying a specific param- 
eter with respect to the baseline parameter set as listed in 
Tableim The baseline values for the sensitivity study cor- 
respond to the default parameter set for ^Snn = 200 GeV 
except for the minimum inter-nucleon separation dis- 
tance which is here set to c?,„i,i = 0.4 fm to match the 
default value in HIJING [4l|. We study the variation of 
all the main Glauber parameters except for the atomic 
mass number. The nuclear radius [R) is varied by ±2%, 
the nuclear skin depth (a) by ±10%; both variations are 
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FIG. 3: Comparison of participant versus binary weighted 
participant eccentricity, as well as the x = 0.13 mixture, as 
a function of A'^part in Au + Au (upper set of curves) and 
Cu + Cu (lower set of curves) collisions at ^Snn = 200 GeV. 
The inset shows the ratio with respect to the participant- 
weighted eccentricity. 



FIG. 4: Comparison of participant eccentricity as a function 
of A^part in Au + Au (upper set of curves) and Cu + Cu (lower 
set of curves) collisions aty^Sp™ = 200 GeV for point-like, 
hard-sphere, and Gaussian matter distributions. The inset 
shows the ratio with respect to the point-like participant ec- 
centricity. 



several times larger than the systematic error assigned 
to their measurements by the authors of Ref. [i^l- The 
nucleon-nucleon inelastic cross-section ((Tnn) is varied by 
more than the experimentally spanned region at RHIC ^ . 
The minimum inter-nucleon separation distance, which is 
not known experimentally, is varied by ±100%. As one 
can see in Fig. [21 both eccentricity definitions are quite 
stable within the studied range of Glauber parameters. 
In particular, this is true for the participant eccentricity. 
Not unexpectedly, the difference between the two defini- 
tions, Epart and erp, is most pronounced for Cu + Cu. 
We have also found that varying two parameters at the 
same time does not increase significantly the observed 
variation in eccentricity. 



B. Particle Production Models 

1. Binary Collisions versus Participants 

The observed particle multiplicity at mid-rapidity 
scales somewhat more strongly than linearly with the 
number of participating nucleons (4^ . This can be 
parametrized by postulating a second (smaller) contribu- 
tion to particle production that scales with the number of 
binary nucleon-nucleon collisions [i^. The MCG model 
as described above, can be extended to implement matter 
distributions produced according to a "two-component" 
scenario, where some of the matter is generated propor- 
tionally to the number of binary collisions. To add this 



feature, it is necessary to define two origins of matter 
production: 

• Participant nucleons, which create the matter by 
means of "exciting" the nucleon. 

• Binary nucleon-nucleon collisions, which create the 
matter locally via a two-body interaction. 

The latter mechanism suggests that the produced mat- 
ter, to be incorporated into the calculation of spatial ec- 
centricity, should be centered between the colliding nu- 
cleons. We achieve this by means of a "pseudo-particle" 
that is located at the center of mass of the pair of collid- 
ing nucleons and keeps track of the "collision-weighted" 
matter with an appropriate weight x. Thus, contribu- 
tions to the eccentricity from participants are weighted 
by while those from binary collisions by x. In 
Fig. [31 the results for purely participant {x = 0) and 
purely collision (a; = 1) weighted participant eccentric- 
ity are shown and compared to the case x = 0.13 which 
has been found to describe the centrality dependence of 
the multiplicity at mid-rapidity according to dA'ch/dr; = 
diVPP/d77 [i^ iVpart +xN,on] El. Independent of cen- 
trality, the collision weighted Cpart values are shifted to 
larger eccentricity for all centralities, similar to what is 
known for the standard eccentricity The presum- 

ably realistic case of a; = 0.13 yields about 10% larger 
eccentricity in the most central Au -|- Au collisions. 



2. Effects of Smeared Matter Distributions 



^ A posteriori, this is justified since the dependence on ctnn turns 
out to be small. Thus, this approach will allow the treatment of 
the systematics at all collision energies in the same way. 



The driving force for hydrodynamic elliptic flow 
is not directly the eccentricity of the distribution 
of participating nucleons or binary nucleon-nucleon 
collisions, but rather the anisotropy of the pressure 
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gradients in the initially produced hot matter after it 
has thermalized. The transverse distribution of this 
matter will be smeared around the transverse positions 
of the participating nucleons or binary collision points. 
To study the behavior of the eccentricity under differ- 
ent definitions, we introduce a general procedure for 
incorporating a variety of matter density distributions. 
In this approach, the contribution of each matter 
production point (i.e. the center of a participant nucleon 
or a binary collision) at (xi,yi,Zi) is smeared according 
to P{x — Xi,y — yi,z — z^), leading to a continuous 
weight function defined at all space-points in the trans- 
verse plane, w{x,y,z) = J2 Pi.^ - Xi,y - yi,z - Zi). 
Averages and higher moments in space-time are then 
calculated (for individual events) using this weight func- 
tion, e.g. {x'} — J X w{x,y, z) dxdydz. The point-like 
MCG cases described above correspond to the choice of 
P{x -Xi,y- T/j, z - Zj) = S{x - Xi)S{y - yi)6{z - Zi). 
We look at two different, azimuthally symmetric, 
parametrizations for the smearing profile 

• Hard-sphere smearing, Phs(?') oc ^(^hs ^ ''^) 



Gaussian smearing, Pg(r) 



exp (-r^/2CT? 



where = (x — Xi)"^ + {v — ViY + (-^ ~ '^iY ^-^id de- 
notes the step function. In the following, we estimate 
meaningful choices for the parameters, i?hs and CTg. For 
i?hs, it makes sense to use the interaction radius, i?hs = 
D/2 = yJcfi^/A'K, since the nucleons in the Glauber 
model are assumed to interact if their centers are within 
the "ball diameter" D, Eq. Q. For 200 GeV, this corre- 
sponds to i?hs ~ 0.6 fm. Then, matching the root-mcan- 
square (RMS) width of the Gaussian distribution to the 
RMS width of the hard-sphere, (t| = R^^l^- leads to 
(Tg « 0.3 fm. In Fig. 2] we show the results for the partic- 
ipant eccentricity calculated for point-like, hard-sphere, 
and Gaussian local matter distributions, using this set of 
parameters. The comparison reveals that for both colli- 
sion systems the way the produced matter is distributed 
around the MCG interaction points does not significantly 
influence the observed value of the eccentricity except for 
extremely small systems. This also shows that the quan- 
tum mechanical uncertainty on the transverse positions 
of the interaction points has no major influence on the 
initial source eccentricity. Note that similarly to what 
was reported in Ref. [43 we find significantly different 
ccntrality dependence for Cpart if we allow smearing out 
the local matter sources to a very large extent. For ex- 



ample. 



-part 



in Au + Au collisions at all centralities does 



not exceed 0.15 for Gaussian smearing with (Tg = 2fm. 



IV. CORRELATIONS AND FLUCTUATIONS 

In this section, we focus on eccentricity cumu- 
lants (which enter the discussion and interpretation of 
elliptic flow data) in the context of fluctuating initial 
conditions and for different realizations of the Glauber 
model initial state. 



A. Sensitivity of the Event-Plane Method 
to Underlying Flow Fluctuations 

There are different ways of extracting the elliptic flow 
from data: the event-plane method, two-particle correla- 
tions, multi-particle cumulants, etc. (see Refs. [TtI. 
Each flow measurement is based on a different moment of 
the final-particle momentum distribution and thus is dif- 
ferently affected by event-by-event flow fluctuations (and 
non-flow correlations). If elliptic flow is proportional to 
the spatial anisotropy, the eccentricity scaling should be 
performed with corresponding moments of the partici- 
pant eccentricity [2^, [20| . 

It has been explicitly stated [2^ (see also [l^) that 
the event-plane method (w2{EP}), used by the PHO- 
BOS experiment to measure elliptic flow, really mea- 
sures (wl) rather than (W2), i.e. the RMS rather than 
the mean of V2- More sp ecific ally, it has been claimed 
that {V2) < {EP} < depending on the event- 

plane resolution, with the upper limit being approxi- 
mately reached under RHIC conditions [43|. Here, (. . .) 
indicates an average over many collision events. In the 
following, we will investigate and confirm this claim. 

In the event-plane method [l3|, one uses the particles 
from one side of the detector (subevent A) to estimate 
the event plane, the plane relative to which the flow de- 
velops, given by ^'^ = ^tan~-^ sin(2(/)i)/ ^ cos(2(/)i)]. 
One then correlates the particles from the other, sym- 
metric, side of the detector (subevent B) with this event 
plane to obtain the uncorrected flow signal for the given 



subevent, 



obs 



{cos {2(j)i — 25*^)}. Here, as before and 
in the appendices, {...} indicates the average over an 
individual event. The roles of subevents A and B can 



then be interchanged to obtain 'I'f 



,,obs 
"2, A' 



and thus the 



observed flow signal for the whole event ^2'^*'. Assuming 
small dynamical and non-flow correlations one has (for 
symmetric subevents) 



i,obs 
^2,B{A) 



{V2) (cos 



2*9 - 2* 



a{b: 



(2) 



where the average is, unlike before, not over particles 
in a given event but over events in the given ccntral- 
ity, r] and pt bin. In this equation, 5*2 stands for 
the actual event plane angle, which deflnes the orien- 
tation of the V2 signal in a particular event [l^l, and 



(cos ^25*2 — 2^!2^^^^ ^ = R quantifies the event-plane 
resolution, which itself depends on V2 The resolu- 
tion can be estimated based on data alone. 



i?= v/ (cos (2^-^-2*3)) (3) 

leading to {vf''') = (ua) R. The presence of a frac- 
tion (/bkg) of uncorrelated background particles in ad- 
dition to the N particles that carry the flow signal corre- 
sponds (restricting ourselves to second harmonic contri- 
butions) to a distribution of 



dN 



(l + /bkg) 



N 
2^ 



2V2 



bkg 



■cos (2(/)- 2^-2) 



, (4) 
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FIG. 5: Dependence of a, Eq. ((6)1, on the event-plane resolu- 
tion, R, for P{v2) uniform with various combinations of 1^2, min 
and "1^2, max 5 number of signal (N) and background (fbksN) 
particles in the MC calculation (see text for more details). 
The shaded band covers the parameter errors obtained from 
a polynomial fit to the data. 



various values of N between 30 and 1000, as well as val- 
ues of /bkg between and 5. For each set of parameters, 
10^ events have been simulated. The number of parti- 
cles in each event is not allowed to fluctuate, i.e. exactly 
N + fh]^gN particles are created in every event. Within a 
reasonable spread that increases with decreasing resolu- 
tion, the a values are found to lie on a common curve as 
a function of R, with no or weak dependence at most 
on the chosen simulation parameters. The PHOBOS 
t;2{EP} measurements lie in the range of 0.15 ^ 0.55 
for Au + Au and 0.13 <i?< 0.33 for Cu + Cu where 
L'2{EP} scales approximately with the RMS of the un- 
derlying V2 distribution. This result is supported by re- 
cent, full detector simulations with HIJING MC events 
that incorporate known, but fluctuating. V2 values. For 
the measured value of (yv2/{v2) ~ 40% [3^ this implies 
that the PHOBOS i;2{EP} measurements are about 10% 
larger than the mean elliptic flow. 

B. Participant-Eccentricity-Scaled Elliptic Flow 



leading to an apparent suppression of the observed flow 
signal by 1 + /bkg- Correcting for this effect, we arrive at 
the final expression for the flow measured via the event- 
plane method 



1'2{EP} 



^'2 



1 + /bkg 

R 



(5) 



This includes the two main experimental corrections, 
namely for suppression and event-plane resolution 

From Eq. ([5]) it is not obvious if the U2{EP} scales with 
the mean or the RMS of an underlying V2 distribution, 
and how such a behavior depends on R. However, one 
can study this behavior numerically with a MC calcu- 
lation that creates events with V2 distributed according 
to a given distribution, P{v2). For every event, we take 
N particles (at mid-rapidity) to carry the flow signal, 
according to ANld(f) oc 1 -|- 2v2 cos(20), and /bkg N parti- 
cles (/bkg > 0) to represent the uncorrelated background, 
which are added with a uniform azimuthal distribution. 
To relate the obtained W2{EP} to a moment of the input 
distribution, we implicitly define the exponent a accord- 
ing to 



V2{EP} 



(6) 



where the ensemble average ((•••)) is calculated with 
the underlying V2 distribution, P{v2). One obtains the 
sensitivity of the event-plane method to P(f2): a — 1 
corresponds to a scaling of W2{EP} with the mean, 
a = 2 to a sca ling with the RMS of P{v2), also often 
denoted V2{2} [4y|. In our calculation, we choose V2 
to be uniformly distributed between 112, min and 'y2,max, 
i.e. the true mean and RMS are 0.5 (t'2,min + V2,ma.x) and 



\/ (^2,min + ^'2,mint'2,max + fi,max)/3, respectively. 

Fig. [5] shows the result of the calculation for various 
combinations of V2 min and W2.max between 0.01 and 0.3, 



It is expected [2l|, [22| that V2/e scales with the trans- 
verse charged particle area density at mid-rapidity ac- 
cording to V2/e oc 1/ S dNch/dy, where S is the overlap 
area in the transverse plane. The data are measured in 
bins of centrality, translating into 



(7) 



|y|<i 



Here, W2{EP} is the PHOBOS estimate of the ensemble- 
averaged elliptic fiow according to Eq. (O, and e{EP} 
is a suitable ensemble-averaged initial source eccentric- 
ity. In view of the discussion of the preceding subsection, 
W2{EP} is about f 2{2}, the RMS of the underl yin g V2 dis- 
tribution. Following the suggestion of Refs. [2J, [23 we 
therefore scale it with e{EP} = epart{2} = (cpart^), the 
RMS of the participant eccentricity distribution obtained 
from the MCG model. The overlap area 



S = TT 



y xy 



corresponds to the area of the tilted overlap ellipse (see 
Appendix 1^1 and especially Eq. (jA8|) ). 

To arrive at the centrality average for epart{2} and S, 
we fold their iVpart dependence with the distribution of 
A^part values obtained for each centrality bin from full 
detector simulations. Independent of species and energy, 
we scale the elliptic flow data by 0.9 and the mid-rapidity 
yields by 1.15 to convert from pseudo-rapidity (77) to ra- 
pidity (y) as described in Ref. [13, Fig. E] shows 
the result for Cu + Cu and Au 4- Au collisions at 62.4 
and 200 GeV. The flow data are from Refs. [1, [Hi, and 
the mid-rapidity yields from Refs. [i^lsH . We distin- 
guish two types of errors that are individually propa- 
gated in the error calculation of the ratios: a) System- 
atic and statistical errors (if available) from data added in 
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FIG. 6; Participant-eccentricity-scaled elliptic flow versus 
transverse charged particle area density at mid-rapidity for 
Cu -|- Cu and Au -I- Au collisions at y^iNN = 62.4 and 
-ysNN = 200 GeV. The horizontal and vertical error bars orig- 
inate from the combined statistical and systematic errors of 
the data (90% C.L.). The shaded boxes are the result of the 
systematic errors assigned to epart{2} and S by the variation 
of the MCG parameters (90% C.L.). 



quadrature to obtain total 90% C.L., and b) systematic 
errors (90% C.L.) assigned to the MC quantities obtained 
by the variation of Glauber parameters with respect to 
the individual baseline values (cp. Table HH). As reported 
earlier [l^jllOl, we find a common scaling between the dif- 
ferent systems. However, within the errors, it is difficult 
to tell whether the almost linear rise of the eccentricity- 
scaled elliptic flow breaks down at larger values of the 
area density which might indicate that the hydrodynamic 
limit is being reached at the top RHIC energy. 



a) 



b) 



c) 




FIG. 7: Schematic of densities for the different approaches: 
a) optical limit, b) full MCG with correlated participants orig- 
inating from each of the two nuclei in one MCG event and 
c) mixed-event MCG with uncorrelated participants where ev- 
ery participant originates from an individual nucleon-nucleon 
collision obtained in a different MCG event. 



C. Cumulants and Correlations 

As mentioned above, it is suggested [1^ that higher 
order cumulant moments of V2 should be proportional to 
analogously defined higher order cumulant moments of 
the eccentricity, including: 

e{2r ^ {e') 

e{Ar EE 2{e'y-{e'). (8) 

In Ref. m Bhalerao and OUitrault (B&O) attempted 
to derive expressions for epart{2} and epart{4} semi- 
analytically, making use of two strong approximations. 
First, the paper contains an implicit assumption that 
all of the participant positions are independent samples 
of some underlying distribution, or at least that any 
correlations between participants do not affect the ec- 
centricity fluctuations [IJ. Second, the expressions in 
Ref. [2^ were obtained using a Taylor expansion, leading 
to a power series in 1/A^part which is then truncated at 
1/iVpait- Based on these approximations, they concluded 



that epart{4} is numerically equal to the standard eccen- 
tricity Eg, vanishing for central (b ~ 0) collisions. This 
would in turn imply that higher order cumulants of the 
flow such as W2{4} are insensitive to fluctuations in the 
participant distribution. In this section we show that 
B&O's assumptions are too strong and that epart{4} for 
Cu -|- Cu collisions differs significantly from eg when bet- 
ter approximations arc made, especially when the role 
of correlations is taken into account, e.g. for the usual 
PHOBOS MCG calculations. 



1. Correlations of Nucleons in the Initial State 

The first approximation in Ref. [23| is to ignore the 
correlations between nucleon participant positions, even 
though we know that there are at least three sources for 
such correlations. 

First, in order to contribute to the produced matter, 
the participating nucleons must hit each other, which 
causes a correlation. For instance, in the case of a pe- 
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FIG. 8: Participant eccentricity and cumulants, tpart, epart{2} and epart{4}, as a function of A^'part in Au + Au (left) and 
Cu + Cu (right) collisions at -^Snn = 200 GeV for the full (upper set of curves) and mixed-event (lower set of curves) MCG 
calculations. The inset shows the ratio of mixed-event to the full MCG results. 



ripheral collision with two or three participants, the over- 
lap region of the nuclei may be something like 3 fm x 1 fm 
but the nucleons will necessarily all be within about 1 fm 
of each other or else they would not be participants. In 
general, the participant positions will tend to be more 
clustered in position space than a random distribution 
since each participating nucleon must hit another one. 

Second, if there are two nucleons from a given nucleus 
which are close together in transverse position, then they 
will have a tendency to be both hit or neither hit. Again 
this will contribute to clustering of participant positions. 

The usual, fuU PHOBOS MCG calculation takes both 
of these effects into account automatically, but the ana- 
lytical expressions, given in Ref. [23| and in Appendix iBl 
here, do not. The two different approximations are illus- 
trated in Fig. [7] (panels a) and b)). 

It should be noted that a precise integral Glauber cal- 
culation at fixed impact parameter should involve a com- 
plete 2{A + i3)-dimensional integral covering all possi- 
ble transverse positions of all nucleons involved in the 
collision [s^, [H, E^l- Such a formulation would in- 
clude the same correlations that are covered automat- 
ically by a MCG. The usual practise of approximating 
the Glauber integral with a single 2-dimensional inte- 
gral (optical Glauber model approximation) is just an ap- 
proximation which neglects all the correlations described 
above. 

A third type of correlation would be genuine nucleon- 
nuclcon position correlations in the wave functions of the 
nuclei that, as mentioned before, generally are ignored in 
all Glauber model calculations. 



S. Uncorrelated Glauber Monte Carlo 

In order to estimate the role of pairwise spatial cor- 
relations among the participant nucleons in a nucleus- 
nucleus collision, we make use of a modified version of 
the MCG code. In the modified version, at first a normal 



event is calculated of which we record the impact param- 
eter and corresponding number of participants. Then, a 
mixed event is constructed from A^part independently cal- 
culated events, which are all required to have the same 
global characteristics, i.e. the same impact parameter and 
number of participants. From each such event, we choose 
one of the participating nucleons, such that in the con- 
structed mixed event none of the participating nucleons 
is correlated to any of the other participants. The result- 
ing approximation of the overlap density is illustrated in 
panel c) of Fig. [71 

Fig.[H]shows the comparison of the full PHOBOS MCG 
calculation with the mixed-event MCG calculation for 
the participant eccentricity and its first two cumulants. 
We find that the contribution of correlations to epart, 
epart{2} and epart{4} is quite important for the smaller 
system, i.e. it is about 20-45% for Cu-|-Cu and about 5- 
10% for Au -I- Au and rather constant over a wide range 
of centrality. Note that the structure seen in Fig. [8] at 
very low A^part values (see also Figs. I2I3I and 2]) is gen- 
uine for the full MCG and not present in the mixed-event 
case. Furthermore, spatial correlations among the par- 
ticipants in the initial state modeled by MCG are found 
to be less important for the reaction plane and standard 
eccentricity definitions (not shown) . For both definitions 
the uncorrelated cases lead to slightly larger eccentrici- 
ties. The deviation between full and mixed-event MCG 
decreases with increasing centrality and is less than a few 
percent for the Au+Au and about 5-15% for the Cu+Cu 
system. 



3. Cumulants in the Extended B&O Approach 

The second approximation in Ref. (23| is to derive ex- 
pressions for Epart {2} and epart{4} using a Taylor expan- 
sion that leads to a power series in 1/iVpart- In Ap- 
pendix [Bl we advance the calculations to higher orders in 
1/iVpait by generalizing Eq. (11) from Ref. [2J|, and show 
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FIG. 9; Comparison of epart{4} from the full MCG (upper set of curves), the mixed-event MCG (lower set of curves) with 
the semi-analytical approach, Eq. (|B37}, and B&O's approximation (Eq. (13) from Ref. [l^), as well as with the standard 
eccentricity as a function of A'part in Au + Au (left) and Cu + Cu (right) collisions at y^s™ ~ 200 GeV. The inset shows the 
ratio with respect to epart{4} from the full MCG calculation. 



how to obtain the analytical terms in a rigorous fashion. 
For £nartl2|, wc obtain the same expression as Eq. (12) 
in Ref. [23, and prove that the O terms are re- 

ally negligible. We also obtain the expansion for epart{4}, 
Eq. (jB37[) . where — in contrast to Ref. — all impor- 
tant terms have been kept. In particular, for central col- 
lisions, when eg — > 0, some terms of O (l/^part) ^^'^ ^'^^ 
negligible and must be kept. 

The values for the ensemble averages over partici- 
pant nucleon distributions (like for example (r^^ or 
(^r''^ cos 20)) in Eq. (jB37p need to be calculated numer- 
ically. We calculate each of these averages as a function 
of ^part using the usual (fuU) PHOBOS MCG code. In- 
serting the numerically evaluated values into Eq. (|B37p . 
leads to the "semi-analytic" result discussed below. 

Fig. [9] shows the results for epart{4}, comparing the 
full PHOBOS and mixed-event MCG with our semi- 
analytical result Eq. (|B37p . and with B&O's semi- 
analytical approximation (Eq. (13) from Ref. [2^ eval- 
uated with the fuU PHOBOS MCG), as weU as with 
the standard eccentricity Cg. For both collision systems, 
our semi-analytical result fully agrees with the mixed- 
event MCG calculation. This is consistent with the fact 
that correlations among the participants are neglected 
in the analytical derivation of Eq. (jB37p . Furthermore, 
it confirms that all numerically important terms have 
been kept in Eq. (|B37|) . The fuU MCG calculation which 
includes participant spatial correlations disagrees with 
the other calculations that neglect them by almost a fac- 
tor of two. Contrary to Ref. [2^, we find that, for the 
Cu+Cu system, epart{4} calculated in the semi-analytical 
approach does not agree with eg, in particular for very 
peripheral and near-central collisions. More important, 
however, is the aforementioned effect of the neglected cor- 
relations. For the Au-I- Au system, epart{4} is found to be 
numerically close to eg (with deviations of less than 10%) 
over a wide range of ccntralities. Only for very peripheral 
and near-central collisions correlations may play an im- 



portant role. For the Cu-|-Cu system, on the other hand, 
epart{4} differs from eg by almost a factor of two over a 
wide range of centralities, implying that correlations can 
not be neglected for the smaller system. 

Fig. (TU] shows a comparison of all eccentricity defini- 
tions used in the present paper (participant eccentricity 
and cumulants, reaction plane and standard eccentric- 
ity), obtained from full MCG calculations with baseline 
parameters listed in Table HI 

As mentioned before, in contrast to the results shown 
in Fig. [51 the authors of Ref. [l^l find in their work that 
epart{4} differs very little from both the standard ec- 
centricity eg and the average reaction plane eccentricity 
(erp), the latter two being almost equal. Recently, the 
authors of Ref. [U have shown that, within a Gaussian 
model of the event-by-event eccentricity fluctuations, the 
identity of epart{4} with (erp) is exact (see Eq. (9) in [5J]) 
as long as the Gaussian widths for erp and for the cor- 
relation term p^y = 2(7^y/(cr^ + "^y) equal. We were 
able to trace the inequality between epart{4} and (eRp) 
present in our MCG model (see Fig. [TO)) to a breakdown 
of the Gaussian model assumptions made in Ref. [U, 
in particular for peripheral collisions and small collision 
systems. We find that in the MCG model the event-by- 
event fluctuations of the correlation term p^y are indeed 
Gaussian for mid-central to central Au -I- Au and most 
central Cu-f-Cu collisions. On the other hand, the event- 
by-event fluctuations of erp arc not well described by 
a Gaussian function for all Au + Au and Cu + Cu col- 
lisions except the most central ones. Furthermore, for 
semiperipheral and peripheral collisions the width of the 
epp distribution does not agree with the width of p^y 
Consequently, for all but the most central collisions, our 
MCG model results are poorly described by the Bessel- 
Gaussian distribution given in Eq. (3) of Ref. [U on 
which the equality of epart{4} and (epp) is based. 
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FIG. 10: Participant eccentricity and cumulants, as well as reaction plane and standard eccentricities as a function of A^part in 
Au + An (left) and Cu + Cu (right) collisions at ^/s^ = 200 GeV obtained from full MCG calculations. 



V. CONCLUSIONS 

The interpretation of the anisotropic flow data mea- 
sured in nucleus-nucleus coUisions at high energy re- 
quires a detailed understanding of the initial source 
anisotropy, which is typically quantified by the eccen- 
tricity of the shape of the nuclear overlap area. In this 
paper, we investigate various ways of defining this effec- 
tive eccentricity using MCG calculations. 

We find that variations in the Glauber parameters 
have only small effects on the participant eccentricity for 
both the Au + Au and Cu + Cu collision systems, while 
the reaction plane eccentricity shows variations on the 
10% level (Fig. [2]). The generalization from participant- 
weighted to collision-weighted interaction point distribu- 
tions leads to an increase in the obtained participant ec- 
centricity, by a constant shift, similarly to what is known 
for the standard eccentricity (Fig. [S]). Over a realistic 
range of parameters, the modeling of smeared matter 
distributions does not lead to significantly different re- 
sults for the participant eccentricity (Fig. [5]). Thus, we 
conclude that reasonable variations in density parame- 
ters, the sources of matter, and their localization have 
only a small effect on the participant eccentricity. These 
results support our initial idea [23 . HOl to use the partic- 
ipant eccentricity definition in conjunction with elliptic 
flow scaling. 

Depending on the event-plane resolution, fluctuations 
in the elliptic flow magnitude influence the measured 
"mean" U2{EP}, and for low resolution bias the mea- 
surement towards the RMS of the elliptic flow distribu- 
tion. For a given event-plane resolution, we find that 
there is a simple connection with the appropriate mo- 
ment of V2 , which appears to be independent of the level 
of uncorrelated background (Fig. [5]). For the resolutions 
achieved with the PHOBOS elliptic flow method [l] this 
gives a scaling with the RMS of V2 for all measured sys- 
tems and energies. We take this into account in the 
eccentricity-scaled elliptic flow (Fig. (5]) by presenting it 



as t'2{EP}/epart{2} together with individual systematic 
errors from data and MC parameter variations. The 
participant-eccentricity-scaled elliptic flow shows an al- 
most linear scal ing with the particle area density as pre- 
dicted in Refs. [2l|, [23l in the low density limit. 



A rigorous attempt to analytically derive non- 
negligible contributions to cumulants of the participant 
eccentricity distribution out to O [l/N^^j.^^ conflrms the 
expressions found in Ref. [l^, for all but epart{4}, where 
our derivation, Eq. jBST}, for the flrst time keeps all lead- 
ing order terms in the series. The numerical evaluation 
of our analytical result for epart{4} agrees with epart{4} 
obtained with the mixed-event MCG calculation, as ex- 
pected, since both ignore all correlations among the par- 
ticipating nucleons. In comparison, the results obtained 
from full PHOBOS MCG calculations, which include spa- 
tial correlations among the participants, imply that pair- 
wise spatial correlations among the participants from the 
collision process itself are quite important, especially for 
the Cu -I- Cu system. For Cu -I- Cu the contribution to 
the participant eccentricity cumulants is about 20-45%, 
while it is about 10% over most of the centrality range for 
Au-I- Au (Fig.[S]). Furthermore, it turns out that epart{4} 
for the Cu + Cu system differs from Cg by about a factor 
of two, while the difference for Au -I- Au is smaller and 
only about 10% for most of the centrality range (Fig. [9l 
Fig. [TU)). Therefore, while correlations among participat- 
ing nucleons may be neglected at the 10% accuracy level 
over a wide range of centralitics in the larger Au-f- Au sys- 
tem, they are crucially important for the smaller Cu + Cu 
system and should not be neglected. These results sug- 
gest that correlations among participants that are cer- 
tainly present in nature and that are to a large extent 
implemented in full MCG calculations are sufficiently im- 
portant that they should be taken into account in any 
study where nuclear geometry is expected to play an im- 
portant role. 
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APPENDIX A: ECCENTRICITY DEFINITIONS 

The spatial anisotropy of the interaction region in the 
transverse plane (in the following as throughout the pa- 
per given by the x- and y-axes) is commonly called "ec- 
centricity" and denoted with the symbol e. It has been 
introduced in Ref. [2l| (called "deformation" , symbol 
5) and in Ref. [l3| (called "spatial asymmetry", symbol 
oix)j while the basic idea of a dimensionless momentum- 
based anisotropy parameter (symbol a) originates from 
Ref. In its most basic formulation, the eccentricity 
definition reads 



The two incoming nuclei, separated by the impact pa- 
rameter 6, can be assumed to be centered at (±6/2,0) in 
the transverse plane such that for a given event the cho- 
sen MC frame coincides with the nuclear reaction frame. 
The "reaction plane" eccentricity is then obtained from 



(A3) 



where = {x^} — {x}^ and cr^ = {y^} — {y}^ are the 
(typically participant-weighted) variances of the nuclcon 
distribution in x- and y-direction in a given event. In 
contrast to eg, which is only defined for the entire en- 
semble of collision events, e^p is defined for each event 
and has its own ensemble average (crp). The reaction 
plane eccentricity can be useful in comparisons with V2 
data where the reaction plane is determined by spectator 
neutrons in a zero-degree calorimeter (23| . Since (crp) is 
numerically very similar to eg (in Au -|- Au the difference 
for A'part > 10 is at most 5%, however in Cu + Cu the 
difference is generally between 15-30%, see also Fig.[TU|. 
it has been used in connection with MCG calculations in 
place of the standard eccentricity 



3. Participant Eccentricity 



^1 



(Al) 



where Bl and R^. characterize the size of the source in the 
X- and y-direction, respectively. Note that this definition 
allows positive and negative values, — 1 < e < 1. It is 
related, but not identical, to the geometrical definition 
of the eccentricity of an ellipse. In this appendix, we will 
summarize the most important definitions of eccentricity. 



The "participant" eccentricity expresses the overlap 
eccentricity in the rotated ( "participant" ) frame (see 
Fig. [T|), denoted by x' and y', which for a given event 
maximizes a'y and minimizes a'^. In principle, the over- 
lap zone will also be shifted with respect to the reaction 
plane frame, but this shift has no impact on the eccen- 
tricity. Generally, the second moments of the position 
distribution in the nuclear reaction plane (or MC) frame 
are described by the covariance matrix. 



1. Standard Eccentricity 



'xy 



(A4) 



Prior to our work [30[ , it has been common practise to 
use smooth, event- averaged initial conditions, for which 
the initial spatial asymmetry in the transverse plane has 
typically been given by the "standard" eccentricity 



(y2) + (a;2) 



(A2) 



where (x^) and (y^) are the second moments of the (typ- 
ically participant-weighted) ensemble-averaged nucleon 
distribution in the x- and y-direction, respectively. We 
follow the notation introduced by Bhalerao & Ollitrault 
in Ref. where (...) denotes an average taken over 
many events (ensemble average), while {• • • } stands for 
an average (over participants) in a single event (sample 
average) . 



where cr^, and a^y ~ {^v} ~ {x}{y} are the per- 
event (co-)variances of the underlying (typically partic- 
ipant weighted) nucleon distribution in the transverse 
plane, given in the original frame. The participant frame 
corresponds to the frame in which E is diagonal. Since 
S is a real symmetric matrix, its diagonalization can be 
accomplished by finding the eigenvalues A that satisfy 
dct(I] — XI) = 0, leading to a second order polynomial in 
A with two solutions: 



al ± 



'^'^ly 



(A5) 



* In Refs. [24l l30l | the reaction plane eccentricity was instead 
named e^td or estandard- 
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These two values of A correspond to a'^ and a'y with the 
larger value (A"*") corresponding to the y' and the smaller 
value (A^) to the x' direction, by definition. This leads 
to the expression for the participant eccentricity [23, [s^l , 



^'2 



^'2 



'^part 



4cr2 



(A6) 



Like eRp, the participant eccentricity is defined on an 
event- wise basis, however in contrast to the previous def- 
initions of the eccentricity, it is non-negative, covering 
the range < Cpart < 1 by construction. The participant 
frame is tilted event-by-event by an angle of ^'part with 
respect to the reaction plane, where 



tan^f 



part 



A 



- (= 



A+ 



r). 



(A7) 



It should also be noted that since the overlap ellipse is 
generally tilted, its area is not proportional to OxOy as 
often assumed, but rather given by 



S 



g'^g'^ — g'^ 



(A8) 



Numerically the ratio, J" J' , is very similar for the Cu + 
Cu and Au -|- Au systems at the same iVpart , larger than 
0.75, and increasing with increasing centrality so that for 
A^part > 20 it is larger than 0.9 in 200 GcV collisions. 



APPENDIX B: CALCULATING CUMULANTS 

In Rcf. m Bhalcrao & OUitrault (B&O) derive the 
behavior of various cumulant moments of the participant 
eccentricity analytically, making use of two questionable 
approximations. First, the paper contains an implicit 
assumption that all of the participant positions are in- 
dependent samples of some underlying distribution, or 
at least that any correlations between participants are 
unimportant to the eccentricity fluctuations [i^. Second, 
the paper uses a Taylor expansion which leads to a power 
scries in 1/A''part which is then truncated at O (1/Apart) 
without proof that the nominally higher order terms are 
actually smaller. 

In this appendix, we will extend the B&O results for 
the case in which any correlations between participant 
positions are still considered negligible, but keeping all 
important terms of the Taylor expansion. In particular, 
we will generalize Eqs. (11)~(14) of Ref. and com- 
ment on Eq. (8). 



1. Generalizing B&O Equation 11 

Following B&O, we will assume that for each event 
N = Apart participants are generated independently 
from an arbitrary underlying 2-dimensional distribution. 



The averaging symbol (/) denotes the average of the 
quantity / over the underlying distribution and/or the 
ensemble average value taken over a large number of 
events. In order to investigate fluctuations, we must 
also consider event- wise averages {f} = jfJ^iLi fi- The 
event- by-event fluctuations are given by i5/ = {/} — (/). 
For convenience of calculation, let us also define / = 
/-(/) such that 5/ = {/}. 
Obviously, 



{Sf) = {f)^0. 



(Bl) 



Next, we evaluate (SfSg) by exhibiting the ensemble 
and event averages explicitly: 



1 



( 



1 



A„,A2 



n— 1 i— 1 



N N^. 



9i,(n) 



\ 



1 



( 



N 



N 



N N 

i=i 

\ 

i—1 n—1 

m J 
\ 



N JVe 



1 = 1 



V 



i=l 



(B2) 



The sums are all finite and their order can be inter- 
changed freely If the participants are numbered ran- 
domly, then (/i) = (/). For example, the average of the 
values for all "participants number 7" over all events 
will just be (x^)- Correspondingly, (/;) — (/) — 0. For 
{fi9i) we have 

im = iifgh) = ifg) ■ (B3) 

For {fiCjj) with i ^ j we have 

(A9,> = (/.>(5,) = (/>(.9>=0, (B4) 

since, following B&O, the positions of participants i and 
j [i 7^ j) in each event are assumed to be uncorrelated. 
It is this last step in Eq. (|B4[) which fails when there 
are correlations between the locations of different par- 
ticipants. Neglecting such correlations, as in B&O, one 
finally arrives at 

^^f^9) - ^iN{fg)+N{N-l){f){g) 



Af2 

M 

N 

ifg) 



N 



if) (g) 



N 



(B5) 



in agreement with Eq. (11) from B&O |23j . 

Generalizing the above derivation to higher orders, we, 
as well as OUitrault [l^l , find that the correct generaliza- 
tion of Eq. (|B5|) for 6^ terms is given by 



fgh 



Af2 



(B6) 



15 



Compared to the S'^ terms, Eq. (|B5|1 . this is suppressed by 
a factor 1/iV. Starting with the (5^ term, the expressions 
begin to get more comphcated. In particular, 



(SfSgShSu) = 



(fghi 



{N-l) 



iV3 ATS 

{fh){9u) + {fu){gh) 



{fg){hu) 



(B7) 



which is actually 0(l/N^), i.e. the same order in the 
number of participants as the term. 

The fifth and sixth order terms can be calculated sim- 
ilarly. Any terms involving single powers like (/) will 
again vanish. The nonzero terms are 



{SfSgShSuSy) 



N - 1 



{fgh) (uv) 
' + {fhu){gv) 



{fgu) (hv) + {fgi') (hu 
Cfhv){gu) + {fuv){gh) + {ghu){fv) 
{ghv){fu) + (guv) if h) + {huv){fg) 

(B8) 



and 

{6f6g6h6uSvSu 



1 

+ {N -l)[{fghu)(vw) 



fghiivw) 



+ {N -l)[{fgh){uvw\ 
+ {N -1){N -2) X 

(/g) {hu) (vw) + 14 perms 



14 perms. 
9 perms 



(B9) 



In general, the dominant terms in the 1/A^ expansion 
should be those composed of products of bilinears like 

l^fg'^ in the case of even powers of d, and those with 

bilinears and one trilinear in the case of odd powers 
of 5. So, we have O (J^"-!) = 0{1/N'') for n > 1 
and ©((52") = 0(l/iV") for n > I. This means that 
we must consider terms up to O (5^) if we want to find 
all O (l/A^^) terms in the series truncated by B&O, and 
terms up to O ((5^) in order to capture the leading be- 
havior in the limit Eg ^ (see below) . 



2. B&O Equation 8 

The Taylor expansion which leads to Eq. (8) in B&O is 
not applicable everywhere as it implicitly assumes that 
l/{Nel) <C 1. Since Es for central collisions, this 
quantity is not guaranteed to be small, and this ex- 
pansion is poorly behaved and formally divergent for 
central collisions. Fortunately, Eq. (8) of B&O is not 
actually needed in order to derive and generalize their 
Eqs. (12)-(14). 



3. B&O Equation 12: Calculating epart{2}^ 



In order to calculate epart{2}^ = (^part)' must first 
express epg^j.^ in terms of the S's. We start with the defi- 
nition, Eq. (|A6|) : 



^part 



Following B&O we have 



(BIO) 

(Bll) 
(B12) 
(B13) 

using {x) = (y) ~ (xy) = {xy'^) = {x'"-y) ~ 0. This leads 
to the exact result: 



{y'}-{yr^{y') + Sy.-S. 

{xy} - {x}{y} = 5.j;y - 5.Jy 



part 



4(5: 



xy 



(r2) 

2es(<5g-<5g) 
(r2) 



(r2)" (r2)" 
25y.^ASl-5l) 



:2\2 



^5151 



(r2 



(B14) 



where (r^) = (a;^) + {y"^)- The second and third terms 
in the denominator are genuinely <C 1 , so it can be safely 
Taylor expanded. The resulting polynomial series is well 
behaved. The leading terms arc and O (cg/A^). All 
terms of O (5^) and higher arc O (eg/iV^) and can be 
dropped since they are at least a full power of 1/7V down 
without any compensating l/eg factors. So, we obtain 



2 26s(5.y2_^2 



2ei6,. 



4(52 

^^xy 



part 



i^r-) {r^} 



(r2 



+ O (eO(53) + O (e°(S4) 



+ 0(e°(55) +0(e[;(56) + ... . (B15) 

This leads to the same result as B&O Eq. (12), except 
that we have also shown that all further terms arc sub- 
dominant: 

1 



(^part) 



2^2 [(l + 3a(r')+4.s(r^cos20)] 



iV2 



(B16) 



Similarly, Eq. (14) in B&O is well-behaved and correctly 
contains all of the leading terms. As noted above, this 
is different for the expansion of Eq. (8) in B&O, which 
does not converge in the limit eg 0. 
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4. B&O Equation 13: Calculating epart{4}'' 

We know from B&O that the O (e, /iV) terms cancel, 
leaving the term as apparently dominant. However, 
in order to confirm this, we need to check that all of the 
nominally higher order terms are actually small. 

In order to organize the calculation, let us write the 
expansion of Equation (|B14p as: 



(B17) 



where A contains all terms of 0{d), B of O {S^), and 
so on. Furthermore, let us define Bq = Iuh^^^qB and 
Co = linic^^oC etc., so that 



.p,,,=So+Co+2?o. (B18) 
Explicitly, we will need the following equations: 



Bo 
Co 



(B22) 



^xy^x ^y 



2(5: 



y^—x-" 



(B23) 



Vn = 



+ ■ 



■^5ly)5l. 



(B24) 



We can now calculate epart{4} : 

= (A') + 2el {B) + 2el (Co) - 2 {AB) 
+ 2el (Vo) - 2 {AC) + 2 {B)^ - (B^) 
+ 4 {Bo) {Co) - 2 {BoCo) + 4 {Bo) (Po) 
- 2 {BoVo) - (C^) (B25) 

where we have kept all terms up to 0{1/N), and the 
leading terms in at O (l/N^) and O (l/N^)- We eval- 
uate the expressions in Eq. (|B25|1 piece by piece, drop- 
ping any terms that would contribute to epart{4}^ at 
0{et/N^) or O (e^'/iV^) for each ?i > 0. Note that 
{r™- cos2n(f>) = ©(e"). This leads to the following ex- 
pressions; 

^ [(l + 3e2)(r4)+4e,(r4cos20)] 



{B) 
(Bo) 



A = 


2€s5y2_^2 2el5r2 




(B19) 




B = 


5l._.,,.+Uly 2e,{5l 


-51) 


4es(5r2(5y2_2.2 


(Co) = 




) 


{r^f 






, '^^{51 + 51) iel5l. 
ir') {r^)' 




(B20) 


{AB) = 


C = 


25y2_^2{5l-5l) 85 


xy^x ^y 


'^^l2_^25,2 


{-Do) - 




85ly5r2 4e,(5,2(52_ 


51)^ 


'iejy2_^2{5l + 5l) 


{AC) = 




ir')' ' {r')' 




ir')' 




6€s5'^25y2_^2 6el5r2 


{51 + t 


1) ^el5% 






ir')' ( 




ir')' 

(B21) 


{B)' ^ 



{BoCo) 
{Bo-Do) 
{CD 



N{r'y 

N{r^y 
1 

N{r'^y 



[{2el + A4){r^) + 8el{r^cos2<j,) 



1 

1 

1 

7V2 



2{r^) 



2es 



^2x3 ,2x3 



2- 



1 



1 

N^{r^)^ 



(r2 

, 2 



(B26) 
(B27) 

(B28) 
(B29) 

(B30) 

(B31) 



10el{rY + 8e,{r^){r^ cos2</)) 

(B32) 

il + 6el){r^)^ + 8es{r^){r^cos2(b) 

(B33) 



{BD = 



1 

7V2 



1 

7V3 

1 

7V3 

1 

7V3 



462-2e2A / +(2 + 146^)^ 



-16e. 



(r^) (r^ cos 20)' 
4(r4)' 8(r4)(r6) 



2{r^) 2(r4)' Q {r^f 



{r^)-2{r^Y 



.2\6 



4(r4) 8{r^y 8 {r'^Y 



(r2 



(r2 



(B34) 
(B35) 
(B36) 
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Assembling, this leads to the final result 



epart{4}' 



+ 



Ar2 



4e, 



(r^ cos 2^) 



16e: 



1 

iV3 



, 2 



16e, 



8(r4)(r6) 



3n 



+ 



o 



O 



ATS 



(B37) 



where we now have all of the leading terms. Terms which 
have been dropped are down from the leading terms by 
at least a full factor of 1/A^ without any compensating 
1/es factor. B&O ^ left out the 0{l/N^) term and 
most importantly the 0{l/N^) term. For Cu + Cu, 
the O [l/N^) term tends to be comparable to the "lead- 
ing" term. For central collisions, as es vanishes, the 
O (l/N^) term becomes dominant and certainly cannot 
be neglected. 



I. Arsene et al. (BRAH MS), Nucl. Phys. A757, 1 (2005), 
|arXiv:nucl-ex/0410020] 

K. Adcox et al. (PHENIX), Nucl. Phys. A757, 184 



(2005) , |arXiv:nucl-ex/0410003 

B. B. Back et al. (PHOBQS), Nucl. Phys. A757, 28 



(2005) , |arXiv:nucl-ex/0410022 
J. Adams et al. (STAR), Nucl. Phys. A757, 102 (2005), 



|arXiv:nucl-ex/0501009 

K. Adc ox et al. (PHENIX), Ph ys. Rev. Lett. 89, 212301 



(2002) , |arXiv:nucl-ex /0204005 
S. S. Adler et al. (PHENIX), Phys. Rev. Lett. 91, 182301 



(2003) , |arXiv:nucl-ex/0305013 

S. S. A dler et al. (PHENIX), Ph ys. Rev. Lett. 94, 232302 
(2005) , |arXiv:nucl-ex/0411040| 

B. B. Back et al. (PHOBQS), Phys. Rev. Lett. 89, 222301 



(2002), arXiv:nucl-ex/0205021 



B. B. Back et al. (PHOBQS), Phys. Rev. C72, 051901 



(2005) , |arXiv:nucl-ex/0407012 
B. B. Back et al. (PHOBQS), Phys. Rev. Lett. 94, 122303 



(2005), arXiv:nucl-ex/0406021 

K. H. Ackermann et al. (STAR), Phys. Rev. Lett. 86, 



402 (2001), [ajXiv:nucl-ex/0009011 

J. Adams et al. (STAR), Phys. Rev. C72, 014904 (2005), 



arXiv:nucl-ex/0409033 



H. Sorge, Phys. R ev. Lett. 78, 
arXiv:nucl-th/9610026| 
H. Sorge, Phys. R ev. Lett. 82, 
arXiv:nucl-th/9812057[ 



2309 (1997), 
2048 (1999), 



P. F. Kolb, J. SoUfrank, and U. Heinz, Phys. Rev. C62, 



054909 (2000), arXiv:hep-ph/0006129 



U. Heinz and P. F. Kolb, Nucl. Phys. A702, 269 (2002), 



arXiv:hep-ph/0111075 

A. M. Poskanzer and S. A. Voloshin, Phys. Rev. C58, 



1671 (1998), arXivmucl-ex/9805001 



J.-Y. Ollitrault, Phys. Rev. D46, 229 (1992). 

P. F. Kolb, P. Huovinen, U. Heinz, and H. Heiselberg, 



Phys. Lett. B500, 232 (2001), |arXiv:hep-ph/0012137 
[20] R. S. Bhalerao, J.-P. Blaizot, N. Borghini, and 
J.-Y. Ollitrault, Phys. Lett. B627, 49 (2005), 



[21 
[22; 
[23; 
[24; 
[25; 

[26; 

[27 

[28; 

[29 
[30 
[31 
[32; 
[33 
[34; 
[35; 
[36 
[37 
[38; 
[39 



arXiv:nucl-th/0508009 



H. Hei selberg and A.-M. Levy , Phys. Rev. C59, 2716 
(1999) , |arXivmucl-th/9812034[ 

S. A. Volo shin and A. M. Poskanz er, Phys. Lett. B474, 
27 (2000), 



arXiv:nucl-th/9906075 



erao and J.-Y. Ollitrault, Phys. Lett. B641, 



R. S. Bha 

260 (2006) 

B. Alv er et al. (PHOBQS), Ph ys 
(2007) , |arXiv:nucl-ex/0610037| 
A. Adare et al. (PHENIX), Phys 



arXiv:nucl-th/0607009 
(PHOBQS) 



Rev. Lett. 98, 242302 



Rev. Lett. 98, 162301 



(2007), [a?Xiv:nucl-ex/0608033 
T. Hirano, M. Isse, Y. Nara, A. Ohnishi, and K. Yoshino, 
Phys. Rev. C72, 041901 (2005), [arXivmucl-th/0506058 
B. Alver et al. (PHOBQS), Phys. Rev. Lett. 96, 212301 



(2006) , |arXiv:nucl-ex/0512016 
G. Roland et al. (PHOBQS) 



(2006) , |arXiv:nucl-ex/0510042 
M. Miller 



and 



R. 



arXiv: nucl-ex/03 12008 



Nucl. Phys. A774, 113 
Snellings (2003), 



S. Manly et al. (PHOBQS), Nucl. Phys. A774, 523 



(2006), arXiv:nucl-ex/0510031 



B. Alver et al. (PHOBQS), J. Phys. G34, S887 (2007), 
nucl-ex/0701054. 

B. Alver et al. (PHOBQS), submitted to Phys. Rev. Lett. 



(2007), arXiv:nucl-ex/0702036 



P. Sorensen (STAR), J. Phys. G34, S897 (2007), nucl- 
ex/0612021. 

M. L. Miller, K. Reygers, S. J. Sanders, and P. Steinberg 



(2007) , |arXiv:nucl-ex/0701025 

A. Bialas, M. Bleszynski, and W. Czyz, Nucl. Phys. 
Bill, 461 (1976). 

A. Bialas, M. Bleszynski, and W. Czyz, Acta Phys. 
Polon. B8, 389 (1977). 

T. W. Ludlam, A. Pfoh, and A. Shor (1986), Brookhaven 
National Laboratory Report - BNL-37196 (Feb. 1986). 
A. Shor and R. S. Longacre, Phys. Lett. B218, 100 
(1989). 

W. M. Yao et al. (Particle Data Group), J. Phys. G33, 
1 (2006). 



18 



[40] 
[41] 
[42] 
[43] 
[44] 

[45] 

[46] 



C. W. De Jager, H. De Vries, and C. De Vries, Atom. 

Data Nucl. Data Tabl. 36, 495 (1987). 

X.-N. Wang and M. Gyulassy, Phys. Rev. D44, 3501 

(1991). 

B. B. Back et al. (PHOBOS), Phys. Rev. C70, 021902 



(2004) , |arXiv:nucl-ex /0405027 
D. Kharzeev and M. Nardi, Phys. Lett. B507, 121 



(2001), [arXiv:nucl-th/0012025 

P. F. Kolb, U. W. Heinz, P. Huovinen, K. J. Eskola, 
and K. Tuominen, Nucl. Phys. A696, 197 (2001), hep- 
ph/0103234. 

W. Broniowski, P. Bozek, and M. Rybczynski (2007), 
larXiv:0706.4266] 

N. Borghini, P. M. D inh, and J.-Y. Olhtrauh , Phys. Rev. 
C64, 054901 (2001), [arXiv:nucl-th/0105040 



[47] 
[48] 

[49] 

[50] 

[51] 

[52] 
[53] 
[54] 



J.-Y. OUitrault, private communication (2006). 

C. Adier et al. (STAR), Phys. Rev. C66, 034904 (2002), 



arXiv:nucl-ex/0206001 

B. B. Back ct al. (PHOBOS) Phys. Rev. C65, 061901 
(2002), [arXiv:nucl-cx/0201005| 

B. B. Back ct al. (PHOBOS), Phys. Rev. C74, 021901 



(2006) , |arXiv:nucl-ex/0509034 
B. Alve r et al. (PHOBOS ), submitted to Phys. Rev. Lett. 

(2007) , larXiv:0709.4008l [nucl-ex] . 

W. Czyz and L. C. Maximon, Annals Phys. 52, 59 (1969). 
A. Bialas and W. Czyz, private communication (2001). 
S. A. Voloshin, A. M. Poskanzer, A. Tang, and G. Wang 
(2007) , larXiv:0708.0"800] 



